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ABSTRACT 

Historically the velocity scatter seen on local Hubble plots has been attributed to the peculiar 
velocities of individual galaxies. Although most galaxies also have uncertainties in their distances, 
when galaxies with accurate distances are used recent studies have found that these supposed 
peculiar velocities may have preferred, or discrete, values. Here we report the interesting result 
that when these discrete components are identified and removed from the radial velocities of the 
SNela galaxies studied in the Hubble Key Project, there is evidence for a residual oscillation, or 
ripple, superimposed on the Hubble fiow. This oscillation has a wavelength near 40 Mpc and, 
because its amplitude is small compared to that of the scatter in velocities, it becomes visible 
only after the discrete components are removed. This result is interesting because even if this 
ripple has been produced by a selection effect, the fact that it is only revealed after the discrete 
velocities are removed implies that the discrete velocities are real. Alternatively, if no selection 
effect can be identified to explain the ripple, then both the discrete velocities and the ripple 
together become very difficult to explain by chance and these results could then have interesting 
cosmological consequences. 
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nents in galaxies introduce a scatter in the Hubble 
plot that is much larger than can be explained by 
the errors of measurement. This scatter has been 
explained historically by peculiar velocities, al- 
though in some cases this has b een questioned sim- 
ply because of their large size (lRussellll2005allbl l3). 
If the quantization found in these velocities is real 
it would appear to rule out the peculiar velocity 
interpretation. 

In this paper we report a remarkable result 
in which a low-level ripple is clearly seen to be 
present in the residual Hubble plot after the dis- 
crete velocities are identified and removed from 
the velocities of the SNela galaxies studied in the 
Hubble Key Project. Because its amplitude is 
small compared to the scatter in velocities it is not 
visible before the discrete velocity components are 
removed. One model that might explain this type 
of oscillation i n the Hubble flow has b een discussed 
previously by iMorikawal (|l99d Il99ll ). 



It has been demonstrated (jTifftl 119961 |1992D, 
(and related papers) that there appear to be dis- 
crete "velocity periods" present in the redshifts of 
galaxies. The most obvious of these was found in 
common spirals and showed discrete velocity com- 
ponents near 36, 72, 145, 290, etc., km s~^. In 
each of the other period groups detected, the ve- 
locities showed this same octave related, or dou- 
bling nature. More recently we have found ev- 
idence that the extrapolation of Tifft's periods 
to higher components, using this doubling rela- 
tion, leads to discrete velocity components that 
appear to be visible in the radial velocities of 
all galaxies whose distan ce s are accurately know n 
()Bell and Comeaull2003allbl : iBell et al.ll2003l \2004 . 



These small non-cosmological redshift compo- 
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2. Discrete Velocities Defined in Previous 
Work 

The discrete "velocity" components in galaxies, 
found by .Tifft (,1996. . , 1997.) and recently confirmed 
by us, can be expressed by the relation: 

ZiG[N,m] = iZiG[N,nmax]/2"')rn=0,1.2,3..oo 

(1) 

Here the integral values of the quantum num- 
ber N correspond to the different "velocity" pe- 
riods identified by Tifft, and the quantum num- 
ber m represents the number of halvings (via 
2"™) below the relevant maximum intrinsic red- 
shift component z^p \N, Umax] in each iV-group 
( Bell and Comeau 2003a, The maximum for 
galaxies is equal to the minimum found previ- 
ously for quasars in each relavent group (|Bell 
[2002a,b,c,d, 2007). 

Until a physical explanation of the discrete 
components becomes known and it is determined if 
these are discrete velocities as suggested by Tifft, 
or preferred redshifts, in which case they might be 
related to the atom itself, or even some systematic 
effect, it is not possible to explain what the N and 
m quantum numbers might be due to. 

The discrete components in galaxies get quite 
small as m increases and because of this cannot be 
resolved at high m- values. Values for the larger 
(low-m) discrete redshift components in galax- 
ies, and their velocity equivalents, obtained us- 
ing several independent galaxy groups containing 
138 galaxies with accurate Tulle y -Fisher distances, 
are Hsted in iBell and Comeaul (|2003bl Table 4) 
for = 1 to 6. For each value of N the rele- 
N, n„ 



value is given by the highlighted 
Belli (|2002dl . Table 2). The period group 



vant ZiQ 
value in 

found by Tifft to be associated with common spiral 
galaxies, corresponds to the lowest, = 1, group 
of discrete velocities. 

Because we looked at galaxies that were more 
distant than those studied by Tifft we were also 
able to conclude that the discrete components are 
superimposed on top of the Hubble flow. 

3. Analysis of Type lA Supernovae Data 

A preliminary an alysis of the SN cIa data was 
reported previously (jBell and Comeaui,2003b ') and 
it followed the analysis used for spiral galaxies 
(jBell and Comeaull2003allbh where a minimum was 



sought in the RMS deviation in source velocities, 
calculated relative to the nearest discrete veloc- 
ity line superimposed on the Hubble flow. This 
analysis technique is e xplained in more detail by 



Bell et all (|2003ll2004l) . 

In Fig 1 the RMS deviations in Vcmb ve- 
locities, relative to the nearest discrete velocity 
line in Fig 2, are plotted vs H p, using the 3 6 



SNela galaxies listed in Tabic 6 of lFreedman et al 
(j200l[ ). The shape of this curve and a demon- 



stration on how the dip at Ho is produced when 
discrete components are present are discussed in 
detail bv lBell et"al] (|2004l ). A clear best-fit feature 
is visible here at Ho = 58 km s~^ Mpc~^, as was 
found for the 138 spiral ga l axies studied prev iously 



()Bell and Co meaul l2003bl: iBell et al.l l2003l). We 



have shown pre viously ( Bell and Comeaul 2003al : 



Bell et al.ll2003l ) that there is nothing in our anal- 



ysis procedure that can always produce an RMS 
dip at the same value of Ho = 58 in random data. 
The fact that a similar dip at Ho — 58 was found in 
the several independent galaxy groups we looked 
at therefore makes this result very difficult to ex- 
plain by chance. The minimum at Ho = 67 in 
Fig 1 is due to the overall shape of the Vcmb vs 
Distance source distribution. It is not produced 
by the presence of intrinsic components and ap- 
pears in all source distributions, even randomly 
generated ones where there are no discrete veloc- 
ity components present. It can be predicted to be 
slightly lower than the value that would be ob- 
tained by fitting a single straight line to the data 
(as is done in the Hubble Key Project which found 
a value of Ho = 72). 

In Fig 2, the Vcmb velocities of the 36 SNela 
galaxies are plotted vs distance. The slope of the 
discrete velocity lines is determined by Ho and 
their discrete velocity values are defined by equa- 
tion 1 above. Here it can be seen that the RMS 
dip at Ho = 58 is obtained when several sources in 
each [N, to] -group fall along the discrete redshift 
lines over an extended range of distances. This is 
particularly obvious in Fig 2 for the [iV, m] = [2,5], 
[5,8] and [5,9] lines, where 23 of the 32 sources be- 
low 300 Mpc fall along these three lines. Because 
of this it is possible to identify which iV-group the 
sources are from. Unlike the other groups of galax- 
ies studied, at least 60 percent of these sources ap- 
pear to be iV = 2 sources, while the remaining 40 
percent are N — 5 sources. At the same time these 
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sources lie at much greater distances, starting near 
60 Mpc, where the previously studied groups leave 
off, and extending to 450 Mpc. The N values may 
relate to the type of galaxy that generates Type 
la SNe in the redshift range studied here. Tifft 
studied closer galaxies and found that N = 1 and 
2 were most common. However, until an explana- 
tion for the discrete components is found this can 
only be speculation. 

It is also interesting to note that in some re- 
spects the plot in Fig 2 is similar to that in Fig 2 
of Bell (2007.) . where intrinsic velocity or redshift 
components appear to increase with more distant 
objects. 



CO 300 

5 

or 




Fig. 1. — RMS deviation in Vcmb about the in 
trinsic redshift grid lines in Fig 2 vs Hp for Type 
lA Su pernovae. Data are from iFreedman et al 
(l200lh 



4. Hubble Plot After Removal of Discrete 
Velocity Components 

Table 1 lists the SNela sources in col 1, with 
their distances and Vcmb velocities in cols 2 and 
3 respectively. Col 4 gives the intrinsic redshift 
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Data are from IFreedman et al.l (|200ll ) 



Fig. 2.- 
galaxies. 

Dashed lines are N — 2 group and dotted lines 
are N = 5 group. 



[iV, m] (and its associated "discrete velocity") ob- 
tained for each source in the best-fit situation for 
Ho = 58. Col 5 lists the Hubble velocity Vh after 
removal of the discrete velocities. 

In Fig 3 the Hubble velocities (Vh) are plot- 
ted vs distance. A linear regression on the data 
(all 36 points) gave a slope of 58.26 (std. err. 
0.33) km s~^ Mpc~^. Because there are only four 
sources between 300 and 450 Mpc, and because 
the distance uncertainties may increase in propor- 
tion to distance, only those 32 sources closer than 
300 Mpc have been considered for further analysis. 
It is important to understand here that the prob- 
lem related to the increasing distance uncertainty 
is not one of fitting to the ripple, whose amplitude 
may also increase with distance, but one of fitting 
to the discrete velocity lines whose separations do 
not increase. When the distance uncertainties ap- 
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Fig. 3. — Plot of Hubble velocity as a function 
of Distance for SNela galaxies after removal of in- 
trinsic redshifts. 



proach the spacing between these lines the result 
becomes ambiguous. 

Although there are many more SNela galax- 
ies with redshifts between z = 0.4 and 0.8 
measured in the Supe rnova Cosmology Project 
( Perlmutter et al1ll999l ). these have not been in- 
cluded because they lie at redshifts where cosmo- 
logical corrections need to be taken into account, 
and these can be confused with the intrinsic red- 
shifts considered here. 

In Fig 3, where the velocities no longer con- 
tain discrete intrinsic components, there is now 
evidence for a low-level modulation superimposed 
on top of the otherwise linear Hubble slope. It has 
been approximated by the sinusoid (solid curve) to 
which all sources can be reasonably well fitted. 

Since there is nothing in our analysis that could 
have produced this modulation, it appears to in- 
dicate that there is a low-level, sinusoidal oscil- 
lation in the expansion rate of the local Universe. 



Fig. 4. — Plot of residual velocity as a function 
of distance for SNela galaxies after removal of in- 
trinsic redshifts and Hubble slope of 57.9 km s~^ 
Mpc~^. Redshift peaks are located at 30, 70, 110, 
150, 190, 230 and 270 Mpc. 



Changes in the expansion rate are not new, since it 
is now well accepted, at least over a much longer 
time scale, that the expansion rate of the Uni- 
verse is currently accelerating. Investigators have 
argued that ripples in the Hubble fl ow might some- 
day be detected (jPavidsonl |2004 see for exam- 
ple) and others have already reported periodic 



densi t y clumping with redshift (jBroadhurst et aL 
Il990l : Imo et al.lfl992l: iHartnett and Hiranol l2008l) 



(see below for further discussion). 

In Fig 4, a slope of 57.9 has been removed from 
the Vh vs distance plot in Fig 3. Two velocities 
have been plotted for SN 1993ae, which is located 
at a distance of 71.8 Mpc. The open circle is the 
velocity obtained if this source is fitted to the near- 
est N ~2 group (N, m — 2,6). For this assignment 
it clearly does not fit the curve. However, if this 
source is fitted to the nearest grid line (iV, m = 
1,4), it fits the curve well. It is therefore assumed 
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Table 1: Parameters for Supernovae Galaxies. 



Supernova 


D (Mpc) 


VcMB (km s-^)[l] 


Transit.(Disc.Vol.)(km s-^)[2] 


Vh (km s-^)[3] 


SN 1990O 


134.7 


9065 


ZiG [2,6] (1448) 


7617 


SN 1990T 


158.9 


12012 


ZiG [2,5] (2890) 


9122 


SN 1990af 


198.6 


15055 


ZiG [5,8] (3537) 


11518 


SN 1991S 


238.9 


16687 


ZiG 2,5 (2890 


13797 


SN 1991U 


117.1 


9801 


ZiG [2,5] (2890) 


6911 


SN 1991ag 


56.0 


4124 


ZiG [5,10] (888) 


3236 


SN 1992J 


183.9 


13707 


ZiG [2,5] (2890) 


10817 


SN 1992P 


121.5 


7880 


ZiG [5,10] (888) 


6992 


SN 1992ae 


274.6 


22426 


ZiG [2,4] (5742) 


16684 


SN 1992ag 


102.1 


7765 


ZiG[5,9](1772.5) 


5992 


SN 1992al 


58.0 


4227 


ZiG [5,10] (888) 


3339 


SN 1992aq 


467.0 


30253 


ZiG 2,5] (2890 


27363 


SN 1992au 


262.2 


18212 


ZiG [2,5] (2890) 


15322 


SN 1992bc 


88.6 


5935 


ZiG 5,10] (888) 


5047 


SN 1992bg 


151.4 


10696 


ZiG [5,9] (1772.5) 


8923 


SN 1992bh 


202.5 


13518 


ZiG 5,9 (1772.5) 


11745 


SN 1992bk 


235.9 


17371 


ZiG [5,8] (3537) 


13834 


SN 1992bl 


176.8 


12871 


ZiG [2,5] (2890) 


9981 


SN 1992bo 


77.9 


5434 


ZiG 5,10] (888 


4546 


SN 1992bp 


309.5 


23646 


ZiG [2,4] (5742) 


17904 


SN 1992br 


391.5 


26318 


ZiG [5,8] (3537) 


22781 


SN 1992bs 


280.1 


18997 


ZiG [2,5] (2890) 


16107 


SN 1993B 


303.4 


21190 


ZiG [5,8] (3537) 


17653 


SN 19930 


236.1 


15567 


ZiG[5,9](1772.5) 


13795 


SN 1993ag 


215.4 


15002 


ZiG [2,5] (2890) 


12112 


SN 1993ah 


119.7 


8604 


ZiG[5,9](1772.5) 


6832 


SN 1993ac 


202.3 


14764 


ZiG [2,5] (2890) 


11874 


SN 1993ae 


71.8 


5424 


ZiG[2,6] (1449) [1,4] (1158) 


3975(4266) 


SN 1994M 


96.7 


7241 


ZiG [5,9] (1772.5) 


5469 


SN 1994Q 


127.8 


8691 


ZiG [2,6] (1449) 


7413 


SN 1994S 


66.8 


4847 


ZiG [5,10] (888) 


3959 


SN 1994T 


149.9 


10715 


ZiG[5,9](1772.5) 


8943 


SN 1995ac 


185.6 


14634 


ZiG [5,8] (3537) 


11097 


SN 1995ak 


82.4 


6673 


ZiG[5,9](1772.5) 


4901 


SN 1996C 


136.0 


9024 


ZiG [2,6] (1449) 


7575 


SN 1996bl 


132.7 


10446 


ZiG [2,5] (2890) 


7556 



[1] Total measured velocity including discrete component. 
[2] Discrete component obtained. 

[3] Hubble velocity after removal of discrete component. 
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Fig. 5. — VcMB plotted versus Distance before 
(open circles) and after (filled circles) discrete 
components are removed. As can be seen, the rip- 
ple amplitude is too small to be detected before 
the intrinsic values are removed. 

that this one source is an iV = 1 group object. 

Fig 4 clearly shows the presence of a sinusoidal 
fluctuation in the residual velocity (Vr) with in- 
creasing cosmological distance. The best-fit wave- 
length of this modulation is Ag = 39.6 Mpc (period 
= - 1.3 X 10^ yrs, frequency Vg - 2.4 x 10"^^ Hz). 
Its best-fit phase for a sine wave with 39.6 Mpc 
wavelength is 3.18 radians, with redshift peaks at 
29.7, 69.7, 109.7, etc., ±0.3 Mpc. We also note 
that the sources used here are located in all direc- 
tions on the sky, confirming the isotropic nature 
of these results. The amplitude of the oscillation 
may also increase with distance, at least initially. 

The observed ripple, with a wavelength near 
40 Mpc, is si milar to one of the pe r iods r eported 
previously by iHartnett and Hiranol (|2008l ). Their 
distance spacing of (31.7±1.8)h~-'^ corresponds to 
44±8 Mpc for Ho = 72±8 (h = G.72±0.08), which 
is the relevant value for Hq to be used if the dis- 
crete components have not been removed. The 



spacing they find between peaks for this compo- 
nent is thus in good agreement with our 40 Mpc 
ripple. 

The ripple cannot be detected before the intrin- 
sic components have been removed, but it can eas- 
ily be seen by eye after the discrete components are 
removed, both before and after the Hubble slope is 
removed (Figs 3 and 4). There is no obvious sign 
of this ripple before the discrete components are 
removed because the scatter in velocities before 
removal of the discrete components is significantly 
larger than it is after removal of the discrete com- 
ponents. This is demonstrated more clearly in Fig 
5 where Vcmb is plotted versus distance, both be- 
fore and after removal of the intrinsic components. 
There is no obvious sign of the ripple before the 
intrinsic components are removed. 

It is important to point out that to our knowl- 
edge no one has previously reported a velocity os- 
cillation. Previous results using much larger sam- 
ples have found density ripples from a periodic 
clumping with redshift (Broadhurst et al. 1990t 
IMo et al.lll992t IHartnett and Hiranoll2008| r When 
a velocity oscillation is found there is every reason 
to suspect that there might also be density clump- 
ing with the same period. This was pointed out 
by Morikawa (1990), who states that an oscilla- 
tion of the Hubble parameter will cause an appar- 
ent density fiuctuation. If there is a density ripple 
present with a similar period in other source sam- 
ples its detection would not have been prevented 
by the presence of discrete velocity components if 
the distances are accurately known. 

The presence of intrinsic redshift components 
in the radial velocities of galaxies (velocity com- 
ponents that lie only on the high velocity side of 
the Hubble flow as here) will result in a Hubble 
constant that is too high if they are not taken into 
account. This has led to a suggested, revised value 
for the Hubble constant of Hp = 58 km s~ ^ Mpc~^ 
(|Bell and Comeaull2003bHBell et al.ll2003[) . This is 
~ 20% lower than the value reported by the Hub- 
ble Key Project, but in good agreement with the 
value found for intermediate ep och galaxies using 



the Sunyaev-Zel'dovich effect (IJones et al.l 12001 



Mason et alll200ll : iReese et al.ll200a l2003l) . 
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5. Do Selection Effects Play a Role? 

If some identifiable selection effect can be found 
that would produce an oscillation in velocity that 
increases in amplitude as a function of distance it 
would be interesting because it would mean that 
since the removal of the discrete components is 
necessary before the ripple becomes visible, the 
discrete components are then unlikely to be ran- 
dom. If no selection effect can be found it also 
makes the ripple very difficult to explain if it is not 
real since it would then require that two random 
components have worked together in just the right 
way to produce these results by chance, which 
would be highly unlikely. 

6. Discussion 

There are several questions that need to be ad- 
dressed before concluding our investigation. Does 
the presence of the ripple in the raw Vcmb data 
produce the dip seen in Fig 1? Does its presence 
affect the dip? Is the ripple somehow produced by 
the reduction procedure? 

The data analysis used here contains three main 
steps as follows: 

1) The best RMS fit of the data to the nearest 
intrinsic line is calculated as the Hubble slope in 
Fig 2 is varied in unit steps from 45 to 75. These 
values are plotted versus Ho where an RMS dip 
will be detected when data points align with in- 
trinsic components (Fig 1). 

2) The discrete components obtained for the 
best-fit Ho = 58 results are subtracted from the 
Vcmb values to obtain the result in Fig 3. This 
involves a simple linear subtraction that should 
not produce a ripple. 

3) The Hubble slope is subtracted from the data 
in Fig 3 to obtain the residual velocity Vr plot- 
ted in Fig 4. Again this involves a simple linear 
subtraction that should not produce a ripple. 

To investigate whether or not the data reduc- 
tion process produced the ripple we repeated the 
reduction process after first removing the ripple 
in Fig 4 from the raw Vcmb values. This time 
no ripple was detected. The ripple was then sub- 
tracted from the data in Fig 4 and the result was 
compared to the result obtained from the second 
processing. The two were found to be identical, 
as expected. This showed that the ripple is only 



found when one is present in the data and we can 
conclude that the ripple is a real feature that was 
not created by our analysis procedure. 

However, it is also important to show that it 
was not the presence of the ripple in the raw data 
that produced the dip seen at Ho — 58 in Fig 1. If 
the ripple is removed from the raw Vcmb data be- 
fore processing is the dip affected? Does it disap- 
pear completely? Answering these questions will 
also give us an idea of how robust our analysis 
is. In order to investigate this we have included 
in Fig 6 the RMS fit obtained after first removing 
the ripple in Fig 4 from the raw Vcmb data. The 
dip is still clearly present and it is therefore obvi- 
ous that the presence or absence of the ripple in 
the raw data does not affect the fitting of the dis- 
crete components. This is most likely because the 
ripple components are small relative to the separa- 
tion of the discrete components. It also indicates 
that the dip at Ho — 58 is produced entirely by the 
presence of the discrete components in the data. 




100 r \ 

40 45 50 55 60 65 70 75 

Ho 

Fig. 6. — Same as Fig 1 except the sinusoid in 
Fig 4 was removed from the Vcmb data before 
processing. 

However, as can be seen in Fig 6 there are slight 
differences between this curve and the curve in Fig 
1. It can be assumed that there are also small pe- 
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Key Project. These are defined by the same rela- 
tion that also defines those identified previously 
by Tifft and by us in several other independent 
galaxy groups. As in our previous work this re- 
sult is obtained for a Hubble constant of Ho = 58 
kms"^ Mpc~^. We show that when these compo- 
nents are removed from the redshifts of the SNela 
galaxies there is evidence for a low-level sinusoidal 
oscillation superimposed on the Hubble flow. It is 
isotropic in nature and has a period of 40 Mpc. We 
have been unable to identify any systematic effect 
that could have produced this observed oscillation 
in the Hubble flow. If one can be identified it will 
mean that the oscillation is not a real oscillation in 
the Hubble flow. However, it will also mean that 
the discrete components must then be real since 
they have to be removed before the selection ef- 
fect is visible and this is unlikely to happen if the 
discrete components are random ones. If no selec- 
tion effect can be identified that can explain the 
ripple, then both the discrete velocities and the 
ripple together become very diflicult to explain by 
chance and these results could then have signifi- 
cant cosmological consequences. 



Fig. 7.— Mean of Fig 1 and Fig 6. 
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culiar velocities and distance uncertainties present 
in the data. If these affect the analysis differently 
with and without the ripple being present they 
may introduce differences that are truly random. 
If so, their effect can be minimized by averaging 
the curves in Fig 1 and Fig 6. This will minimize 
the effect of random fluctuations but should not 
affect any real features. This has been done in 
Fig 7 where it is seen that while the dip at Ho = 
58 is unaffected, some of the weaker features are 
reduced. We conclude that the dip at Hq = 58 
is a real feature that is likely to have been pro- 
duced by the presence of discrete components in 
the VcMB data as we have argued previously. This 
conclusion is strengthened by the fact that a dip 
has been seen in other source samples at the same 
Ho = 58 Hubble slope. 

7. Conclusions 

We have identified what we believe are dis- 
crete components in the radial velocities of the 36 
Type la Supernovae galaxies studied in the Hubble 
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